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\ Abstract 

| The renormalization of the scalar diquark operator and its anomalous dimension is calculated at two-loop order 

in QCD, enabling higher-order QCD studies of diquarks. As an application of our result, the two-loop diquark 
' anomalous dimension in the MS scheme is used to study the QCD renormalization scale dependence of diquark 

£Nj , matrix elements of the AS = 1 effective weak Hamiltonian. 

^ ' 1 Introduction 

Oh' Four-quark (or tetraquark) qqqq states explain the inverted mass hierarchy of the scalar mesons compared to a qq 
nonet in a variety of theoretical approaches [TJ [3J |31 IH [5] . With the inclusion of a gluonium (glueball) state [5] , the 
i ' scalar spectrum below 2 GeV is then understood as mixtures of gluonium, the qq nonet, and the qqqq nonet. The 
X(3872) 7 and F(4260) [8] mesons can also be interpreted as four-quark states [9]. 

Diquark (qq) clusters are relevant to the internal structure of hadrons (see e.g., [TUl Hi] ) . In particular, Ref. [9] 
uses constituent models for diquark clusters to study four-quark states. The constituent (scalar) diquark masses 
| that emerge in Ref. [5] are in good agreement with QCD sum-rule analyses of diquarks |12[ 113]. providing QCD 
Q\ . corroboration for the diquark model of four-quark states. 

In this paper, we study the renormalization of scalar diquark operators to two-loop order in QCD and thereby 
obtain the two-loop anomalous dimension of the scalar diquark current. As discussed below, the renormalization 
t-H 1 of the diquark operator is an essential component of QCD sum-rule analyses, and the anomalous dimension is also 
necessary for determining the scale dependence of matrix elements of the effective weak Hamiltonian for non-leptonic 
strange particle decays [14]. Our two-loop results thus enable future QCD studies of diquarks to higher loop order. 
£> . The scalar diquark operator in an anti-triplet colour configuration (the "good" diquark in the terminology of 

Ref. 11 j) is given by [12] 



X 

5-H 



J 7 = e Q/ 3 7 Q" (Cj 5 ) l3 q p = e a ^Q a Cj 5 qp , (1) 

where the greek and latin indices respectively represent colour and spin degrees of freedom for the quark fields Q 
and q, and C is the charge conjugation operator. The presence of a transposed quark field in flTJ implies that the 
Feynman rule for the three-point function of the diquark operator and Q, q fields shown in Fig. [1] 

rf = -e a/37 C 75 , (2) 
implicitly transposes the external propagator associated with the Q field. 



Figure 1: Feynman diagram for the tree-level vertex of the diquark operator with the quark fields 
Q and q. The double line represents the Q field that is transposed and the diquark operator is 
denoted by ®. This and all subsequent Feynman diagrams were drawn with JaxoDraw [15] - 
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2 One-Loop Renormalization 



Although the diquark operator is gauge dependent, the theory of composite-operator renormalization |16j implies 
that the diquark operator is multiplicatively renormalizable because there are no lower-dimension operators with 
the same quantum numbers as fT]) LU The one- loop renormalization of the diquark operator can thus determined by 
Fig. [21 which results in the following one-particle irreducible (1PI) Green function for a zero-momentum insertion of 
J 7 in D-dimensions (dimensional regularization) 
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where v is the renormalization scale, the quark mass has been ignored because dimensional regularization is a mass- 
independent scheme, a s = g 2 /{Ait), colour indices have been explicitly shown for the Gell-Mann matrices A a , and a 
covariant gauge with gauge parameter £ has been used. Working in normal (or naive) dimensional regularization^ 
where {7^,75} = [18] in D = 4 + 2e dimensions, and using the (^dimensional) properties of the charge conjugation 
operator CC = —1 and C (7^) C = 7 M [TH] we find 



r (i) 8 r n ,.5 2 1 f d D 



D ( P + kf ( P + kf 







k^kp 
~k?~ 



(4) 



By comparison with the one-loop process determining the renormalization of the scalar current J s — Qq, we see 
that (@| can be related to the (one-loop) 1PI result for the scalar current apart from a numerical factor Cd 
representing the ratio of the different colour factors that occur in the two processes 



r (i) _ I r (o) r (i) =f7 r (0) r (i) 
as represented diagrammatically in Fig. [3] 



(5) 




Figure 2: One-loop Feynman diagram for the renormalization of J 7 . As in Fig. [T] the double line 
represents the (transposed) Q field and the diquark operator is denoted by ®. 




Figure 3: Diagrammatic representation of the relationship ([5]) between two-point functions with 
scalar and diquark operator insertions. The scalar operator is denoted by the solid circle. 

The renormalized diquark operator [J-y] R is defined via the renormalization constant Zd, 
Similarly, the well-known renormalization of the scalar operator is 



\Js 



R 



Zrn J s 



(6) 



(7) 



1 We are grateful for discussions with John Dixon clarifying this point. 

2 We have chosen to work in normal dimensional regularization (as opposed to, e.g., the 't Hooft-Veltman scheme I17H because QCD 
sum-rule analyses of diquarks 1 1 21 1141 have used the normal dimensional regularization scheme. 
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where Z m is the quark mass renormalization constant. Using ([5]) it is easy to see that to one-loop order in the 
minimal-subtraction (MS) and associated schemes 

Zd = z If z H 2 , (8) 

where Z 2 f is the renormalization constant for the quark fields. Landau gauge (£ = 0) is of particular interest in the 
QCD sum-rule analysis of diquark currents, because the Schwinger string used for a gauge-invariant formulation of 
the two-point diquark correlation function vanishes in this gauge |12j . Combining the one- loop Landau-gauge result 
Zip — 1 with ((H)) leads to the one-loop Landau gauge MS-scheme result 

Z d = zW = l + ±°±, (9) 
z tt e 

where we use the dimensional regularization convention D = 4+2e. Eq. (J9j) agrees with the (one-loop) renormalization 
and renormalization-group improvement implicitly implemented in Refs. [121 114] . 



3 Two-Loop Renormalization 

The two-loop diagrams for the renormalization of the diquark operator are shown in Fig. 2) As in the one-loop 
analysis and shown in Fig. [31 each diagram is given by a colour factor Cd multiplying the bare diquark vertex and 
the equivalent diagram with a scalar current. The divergent parts for each of the two-loop diagrams in Fig. [4] are 
expressed in Tabled] in terms of the corresponding scalar diagram T K S [ in the modified minimal-subtraction (MS) 
scheme 



c 



B, 



ie{i,2, ... 11} , 



(10) 



where n/ is the number of active quark flavours and on, and £b are the bare coupling and gauge parameter. A 
number of the Feynman diagrams are clearly related by the exchange of Q and q fields, and hence Table Q] exhibits 
anticipated symmetries r 4 — r 6 , T 7 — T s and T 9 = r u . Note that the colour factors Cd that relate the scalar and 
diquark diagrams are not universally equal to the one-loop result C'd — 1/2, implying that one cannot expect the 
simple pattern of the one-loop result © to persist at two-loop order. The diagrams that are the exception to the 
one-loop pattern (T5 and Tio) require multiple applications of colour algebra identities unique to the Feynman rule 
([2]) ; all other diagrams contain a single application of these identities combined with standard colour algebra factors 
occuring in the renormalization of the scalar operator!! 

The two-loop renormalization procedure first involves the replacement of at and £b with their (one-loop) renor- 
malized expressions (see, e.g., Ref. [21\ ) 



Z a = l + - 

TT 



Z. = l + - 

TT 



33 - 2n f 

4n/ -39 + 9^ 
24e 



a b = Z a a ; 



in the two-loop 1PI Green function 



- r (0) 

1 d — 1 d 



.(2) 



- d ~ x d ■ 

For consistency at two-loop level, (fT3)) requires inclusion of the finite parts of the one-loop calculation §5§ 



3 + 6 



+ 2 (2 + &) - L (3 + &) 



, L = log 


r 2 "i 
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The renormalization constant Zd is then constrained by the requirement that it cancel the divergences in 



Z d Z; 



IF 



,(0) 



,(2) 



where the two-loop MS quark field renormalization constant is [22] 

"£(27+17£) 201 - 12n f + 72f + 9£ 2 



^2F 



l + -f + - 

TT ie V TT 



144e 2 



288e 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 



3 In the previous version of this paper the Table [T] colour factor for diagram 10 in Fig. [4]was erroneous I20| , 
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Figure 4: Two-loop diagrams for the renormalization of the diquark operator where Fi denotes a 
quark loop, T2 a ghost loop and F3 a gluon loop. Implicitly, the Q (double) line extends to the 
insertion of the diquark operator. 



As a benchmark to ensure accuracy in our calculations in Table [TJ we have verified that our results for the scalar 
diagrams lead to the required two-loop MS result Z s = Z m [23] 



Z m — 1 



Q 
7TC 



15 n 



_J_ 
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1 /101 
e Vl8~ 



5ny 
72 



The final QCD result for the two-loop MS diquark renormalization constant is 



1 



a 

7T 
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a^ 2 

7T 
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2880 



64 / 



288 



96 288 



(17) 



(18) 



The cancellation of the L/e terms in Zd that are generated by (jT5J) provides another consistency check on our 

1/2 

calculation. Note that the two-loop Landau gauge result does not uphold the one-loop (£ = 0) pattern Zd = Zm . 
The anomalous dimension for the diquark operator defined by 



v dZd 



Id 



Zd dv 

is easily extracted from (jTSJ to obtain the two-loop MS QCD anomalous dimension for the diquark operator 

2 



a 



a 



ld{a) = 7i- +72 - 



£ 1545 -40n/ 

71 ~ 1 ~ 3 ' 72 ~ 720 



2 16 



In the extraction of the anomalous dimension we have verified that the two-loop coefficients of Zd 

Zd,l , Zd,2 , 



Z,, 



1 



e e 

satisfy the renormalization-group constraint 
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dZ d - 
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(20) 
(21) 

(22) 

(23) 
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Table 1: Results for the two-loop diagrams in Fig. [4] The quantity L = log(— p 2 /v 2 ) and the 
notations for Ai and Bi are defined in Eq. (|10[l , 

where we are working in the conventions of |21) with the (one-loop) /? function and anomalous dimension 5 of the 
gauge parameter given by 



7T 4 O 



(24) 
(25) 



Confirmation of this renormalization-group constraint provides another verification of the accuracy of our results 
given in Table [U 

4 Application and Conclusions 

It has previously been noted that at leading-order, the renormalization scale dependence cancels between the QCD 
perturbative contributions to the diquark decay constants and the AS = 1 effective weak Hamiltonian, although there 
remains some residual scale dependence from non-perturbative terms |14j . As an application of our two- loop results, 
we can explore this scale dependence at next-to-leading order. Following Ref. [14], we consider the combination 



(26) 



where c-(p) represents the renormalization scale dependence of the Wilson coefficient in the AS* = 1 effective weak 
Hamiltonian [24] and <7+(/i) is the scale-dependent scalar diquark decay constant emerging from QCD sum-rules |14j . 
The renormalization-group (RG) factor arising from c_ is [24 



(fi) - exp 



7- (a) da 
/3(a) a 



(27) 
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where in the normal dimensional regularization scheme with rif = 3, the anomalous dimension 7_ (a) i^f| 



7- (") = 
7i = -2 
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7T 
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(28) 
(29) 



Similarly, the anomalous dimension for the diquark operator leads to the following RG factor for the (scalar) diquark 
decay constants 



9+(v>)g+(iJ>) ~ ex P 



-2 



7d(q) da 
/3(a) a 



(30) 



As mentioned above, QCD sum-rule calculations with diquark currents extract gauge-invariant information from the 
two-point correlation function through the insertion of a Schwinger string, which becomes trivial for a line geometry 
in Landau gauge [T2]. Thus for applications to RG behaviour of the diquark decay constants, we use (|2"Tj) with nj — 3 
and £ = 0: 



95 
48 



7i = 1 ; 72 
The resulting RG behaviour of (|26)l is 
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Thus the leading-order cancellation of scale dependence in (|26[) for the perturbative contributions to g + does not 
persist to second order. However, the residual scale dependence associated with (|3"2"|) . which decreases with increasing 
a(/i), does have the right qualitative behaviour to counter the residual scale dependence encountered in Ref. [TJ. A 
more detailed analysis of the residual scale dependence is beyond the scope of this paper because it would require a 
full next-order sum-rule analysis of the diquark decay constants. 

In conclusion, we have determined the MS renormalization constant and associated anomalous dimension for the 
scalar diquark operator at two-loop order in QCD in an arbitrary covariant gauge for normal dimensional regular- 
ization. This result enables future QCD sum-rule studies of diquarks to higher-orders in perturbation theory. For 
example, the divergent terms in the diquark renormalization constant (fl~8| combined with lower- loop O (e) and O (e 2 ) 
terms generate finite parts corresponding to renormalization-induced physical contributions to the diquark correlation 
function. Furthermore, the anomalous dimension of the diquark operator appearing in the renormalization-group 
equation governing scale dependence of the diquark correlation function is an essential feature of QCD Laplace sum- 
rule analyses [22]. Given the relative size of the one- and two- loop terms in (|18l) and pip , these renormalization- 
induced and anomalous dimension effects could be significant in higher- loop extensions of |14) . 
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